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Introduction. Recently, there has been growing interest in the heat and work of nonequilibrium quantum processes [il¬ 
ls^. Studies focusing on this issue have mainly been motivated by an interest in extending the important classical 
fluctuation relations into the quantum regime, e.g., the celebrated Bochkov-Kuzovlev equality [ij and the Jarzynski 
equality (JE)[^. 

Compared with their classical counterparts 32|, which are physically intuitive, the definitions of fluctuating heat 
and work become very delicate in the quantum case. In order to formulate a quantum JE, in a closed quantum system, 
a two-energy measurements (TEM) scheme was proposed by Kurchan 0 to define the work. Although this definition 
still faces criticisms related to the fact that the scheme may destroy the initial quantum-coherent superposition of the 
system , it has been widely accepted in the field [ 13 , [ 13 , [ 13 | . Because closed quantum systems are not common in 
reality, there have also been many attempts to generalize this definition to include open quantum systems [s^ . These 
efforts can be roughly divided into two types of method. The first type [llj, [l3, involves the combination of the 
system and its surrounding heat reservoir as a composite system. The TEM scheme is then conducted on the system 
and reservoir. As the interaction between these is weak, the energy eigenvalue change obtained using the TEM for 
the system is referred to as the internal energy change, while the energy eigenvalue change obtained by the TEM 
for the reservoir is referred to as the heat released from the system, Qtem- Therefore, the work done on the open 
system, Wtem, is the sum of the internal energy change and the heat. The second type of method is based on the 
quantum jump trajectory (QJT) that is unraveled from the Lindblad quantum master equations s i Ea m, 


|24| . I36j| . Under this notion, the energy change of the heat reservoir is continuously measured. This is interpreted 
as the released heat along a trajectory, Qqjt- If one preserves the internal energy change of the system mentioned 
above, an alternative work, Wqjt, is then the sum of the heat and internal energy change along the same trajectory. 
Figure dm schematically illustrates the difference between these two types of work in a two-level quantum system. 

Because very different measurement schemes are involved in the above definitions of work, the nature of the 
relationship between these is not immediately clear. An analogous question related to heat has been raised by De 
Roeck et al. previously. They proved that, if the open system can be described by the Markovian master 

equation, these two types of heat are indeed equivalent. Garrahan and Lesanovsky [33| further emphasized this 
equivalence from the viewpoint of the evolution equation of the characteristic function (CF) of the heat (s^. To my 
knowledge, however, even when using the Markovian setup, the equivalence of the two types of work has not yet been 
explicitly discussed, although intuitively this should be true. A possible reason for this is that, under the action of 
a time- var y in g force, there are no generally valid Markovian master equations except that the driving force is weak 
enough [33, [111 . Very recently, for this specific type of master equation, I developed a CF method [ll] for calculating 
the work Wqjt- It is of note that shortly after the publication of my paper, Silaev et al. [11] proposed another CF 
method for use in calculating the work Wtem for the same quantum system. Hence, now it is important to check 
the equivalence of these two types of work. In this Article, I explicitly show the result. 


Overview of the two CF methods. Let us suppose that the Hamiltonian of a bare system is Hq. This weakly interacts 
with the surrounding heat reservoir, with Hamiltonian Hr, and the inverse temperature /3. The interaction term is 
assumed to be V=S®R, where S and R represent the operators of the system and reservoir parts, respectively. The 
form of V is not the most general form possible, although it is adequate for illustrating our results. Initially, the 
system is in the thermal state, po=exp{—/3Ho)/Ti[e'xp{—/3Ho)]. Then, a weakly driving field is applied from time t=0 
to a final time, </. If one assumes that several approximations are appropriate during the process, which includes 
the weak-coupling, Born-Markov, and rotation wave approximations (RWA), the evolution equation of the reduced 
density matrix of the system p{t) is (setting ?i=l) 

dtpit) = Ctp{t) = -i[Ho + Hi{t),p{t)\ +D[p{t)\, (1) 

where Hi{t) is the interaction term between the system and the field. The H-term denotes dissipation due to weak 
coupling between the system and the reservoir: 


U) 

The operator, S, can be decomposed into a sum of the eigen-operators of Hq, i.e., 

S = Y,S{u:) = Y,S\u:), 

CO CO 


Siuj)ASHu;)-^{S^iu;)Siuj),A} 


( 2 ) 

( 3 ) 


\Ho, S(u)]=-ojS(u}), and S{—ijj)=S'^{uj). These sums are extended over all energy differences, w, of the eigenvalues of 
Hq [331 ■ The rate 7 ( 0 ;) satisfies the detailed balance condition: 


7(-w) = 7(w) exp(-/3a;). 


( 4 ) 
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This condition plays a key role in the validity of the work equalities [1^ and in the following discussions. Master 
equation © is widely utilized in quantum optics, e.g., in describing resonance fluorescence [35l| . 

Equation o can be unraveled into the QJT form for the state vector 1^.1^ EH. By applying this notion, one 
may intuitively define the work Wqjt along each individual trajectory [ H , as described in Fig.(IT|)(b). I have 

shown that the probability density function (PDF) of the work can be calculated using its CF 

^QJT{^^) = = Tro[ifo(0;M)po], (5) 

where /i is a real number, and the symbol E[- ■ ■ ] denotes an average over all QJTs. The operator Ko{t'; u) {0 < t' < tf) 
therein satisfies 




where Ko{tf; /r)=/ is the identity operator, C^, is the dual superoperator of Ct', and 


D*[A] = ^ 7(w) S\uj)AS{u:) - i {S\uj)S{u:), A} 


( 6 ) 

(7) 


Equation ([5]) corresponds to the backward time, t'. Hence, I refer to this as the backward equation. This is a terminal 
value problem rather than a common initial value problem. On the other hand, for use with the same master equation, 
Silaev et al. [1^ presented another CF for use with the work, Wtem- 

^TEM (m) = Tro p{tf;fi)], (8) 

where a modified reduced density matrix p{tf;ij) that satisfied the following expression [l^ was introduced: 


dtp{t;p) = Ct{p)p{t;p) 

= -i[Ho + Hi{t),p{t;p)] 




E>^-Sic.)p{t;p)SHu;) - - {S\uj)Sic.), pit; p)} 


(9) 


with an initial condition of exp[—In order to distinguish this equation from Eq. (jb)), I refer to this as the 
forward equation, since it relates to the forward time, t. The forms of Eqs. ([S]) and (|S]), and Eqs. dH]) and ([5]) appear 
to be very distinct. 


Directly proving the equivalence of the two methods. First, I introduce an alternative operator, 

Ko{s;p) = eKoit';p)E'^^°Q\ ( 10 ) 

for which 0 is the time-reversal operator and a parameter s=tf — t'. By substituting this expression into Eq. ([51), I 
convert this equation into an initial value problem: 


dsKois;p) = Csih)Kois; p) 

= -i[Ho + Hi{s), Kois; p)] 

+ \e-^^^SHuj)Kois;p)Siuj) - i 


( 11 ) 


where i/i(s)=0iFi(t/ — s)0l, and the initial condition iFo(0; p) is equal to exTp[—ipHo]. I have used the detail balance 
condition and assumed that the eigenoperators, S'(a;), are time-reversible. The formal solution to Eq. (Hill is given by 
the following expression: 

Kois;p) = Gis,0;p)[Koi0;p)], (12) 


where (5(5,0;/r)=T_ exp J^jCs'{p)ds' 


is a superpropagator and r_ denotes the chronological time-ordering operator. 


The superoperator Csip) possesses an important property: 


Csifi)A = et,,ip)i&^Ae)e\ 


(13) 
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0: 1st measurement 


t^: 2nd measurement 








W -s -s +coN -oN 

QJT n m + - ^ ^ 



FIG. 1. A schematic diagram describing the two definitions of work based on the TEM scheme on (a) the combined two-level 
qnantnm system (the blue circles) and reservoir (red squares) and on (b) the QJT. In (a), the Hamiltonian, Ho, of a bare system 
has a discrete eigenvector and eigenvalue: Ho|£„)=£„|e„). The Hamiltonian Hr of the reservoir has the following eigenvector 
and eigenvalues: Hr|Cr )=Cr ICr }■ The evolution of the wave vector, of the composite system is unitary under the whole 

Hamiltonian. The green arrow on the right-hand side denotes the projected energy measurements on the system and reservoir 
at time, tf. In (b), the dashed lines denote the Schrodinger-like evolution of the state vector, '(/>o(t), of the system under a 
non-Hermitian Hamiltonian, while thin lines represent occasional jumps of the state vector due to emitting or absorbing a 
quantum lj of energy into or from the reservoir. These energies are recorded by the detectors labeled by the letter D [^. I 
assume their total numbers to be N+ and N-, respectively. Notice that here the projection at time tf is conducted solely on 
the bare system; see the green line. 


where (/i) is the dual superoperator of in Eq. ([Q]): 

CUf,)A = i[Ho + ^ liu:) SHuj)AS{u;) - J {s\io)S{u;), A} 


(14) 


The validity of Eq. (IT^ can be easily proved. Using these operators and the superoperator property mentioned above, 
I rewrite Eq. using the following expressions: 


= Tr[0tHo(t/;Ai)0e-*'^"Vo] 

= Tr[0tG(t/, 0; Ai)[iTo(0; Ai)]0e-*'^^Vo] 
= Tr[G*(0, t/;/x)[0tHo(O;/r)0]e-*'^^Vo] 
= Tr[e*^^'>G(t/,0;M)[e-*'^^Vo]], 


(15) 

(16) 

(17) 

(18) 


, and T+ denotes the antichronological 


where G(t/, 0; fi) is the dual superoperator of G*(0, t/; fJ,)=T+ exp C*{fi)dT 
time-ordering operator. It is clear that G(t/,0;/x) is in fact the superpropagator contained in Eq. Therefore, 
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these two CF methods are equivalent, i.e., ^QJT{^J)=^TEM{^J)■ Since the backward and forward time parameters 
are involved, the current situation is very similar to the case of the Kolmogorov backward and forward equations 
employed in classical stochastic theory [^, . 


Microscopic basis for the backward equation. Equation has a microscopic origin. This expression was obtained 
by reducing an equation relating to the CF of the work, Wtem, defined for the composite system [H, [H, into 
the degrees of freedom of the system. The above proof implies that Eq. shall be derived in an analogous manner, 
even though this expression was obtained solely by employing the notion of QJT . After all, the quantum master 
equation and its manner of unraveling can be thought of as effective theories. However, if one follows the forward 
approach, as employed by Silaev et ah, it is impossible to directly arrive at the backward equation ([6]). Let us recall 
another CF method for computing the PDF of the work, Wtem, of the composite system. First, I explicitly write its 
whole Hamiltonian, H{t')=Ho+Hi{t')+Hr+V. There exists the following evolution equation relating to the operator, 
K{t'-p) HI: 

dt’K{t'-p) = -i[H{t'), K{t'-p)] - K[t'-p)i[H{t'), (19) 

with a terminal condition of Ar(0; p)=I. The CF, $t£;m(m), equals Tr[K(0; p)po®Pr\, where pr=exp(—/3iLr)/Tr[exp(—/SiL^)] 
is the canonical density matrix of the reservoir. A brief explanation of Eg. (1191) is given in Appendix I. This equation 
appears to be very similar to Eq. ([6]). In particular, as $T£;M(M)=Tro[Trj.[Ar(0;/r)pr]po]j one may naturally assert 
that the latter could correspond to the reduced effective equation of the former, while the term Tiv[- • •] (the trace 
over the reservoir) could correspond to the previous Ko{0-,p). In the following discussion, these two conjectures are 
verified. 

Let us introduce the time evolution operator Uor{t')=exp[—i{Ho + Hr)t'] and rewrite Eq. (ITOl) using the interaction 
picture, 

dTKiit'; p) + i[Hi{t'), Ki{t'-p)] + iKiit'-, p)e^^^^o , H/(t')] 

= -i[Vi{t')Kiit'; p) - Ki{t'; p)Vi^{t')], (20) 

where the operator subscripts / denote that these are the interaction picture operators, and 

Vif,(t') = ( 21 ) 

Note that I have moved all terms that do not involve the interaction term, V, to the left-hand side (LHS) of the 
equation. This is in preparation for the perturbation calculation below. Eq. (1201) has an integral form: 


Ki{t']p) = I+i / dT[Vi{T)Kj{T-,p) - Kj{t;p)Vi^{t)] 


ftf 


ftf 


iU-Ho 


+i J dT[Hl{T),Ki{T; p)]+i J dTKi(T\p)e'^ °[e 

I substitute Eq. (1^ into the right-hand side (RHS) of Eq. (OIH) and obtain 

dt-Ki{t'; p) + i[Hiit'), Ki{t'- p)] + iKi{t '-, Hi{t')] 
= r dTViit')[Vi{T)KiiT;p) - KiiT;p)Vi^iT)] 




( 22 ) 


/ dT[ViiT)Ki{T;p) - KiiT;p)Vi^iT)]Vi^it') - f[IA(t') - lA^(t')] 
Jt' 

dTVi{t')[H[{T),Ki{T;p)]- dT[H({T),Ki{T;p)]Vifj,{t') 

f ' dTViit')Kj{T-p)e^>^^°[e-^>^^°,H[{T)] 

Jt' 

r dTKj{T- p)e^>^^° , HliT)]Vj^{t’). 

Jt' 


(23) 


Multiplying both sides by pr and taking a trace over the reservoir, I transform the LHS of the above equation into 
LHS = dt'K^i{t'-p)+i[Hi{t'),K^i{t'-p)]+iKoi{t'-p)e^^^'^[e-^^^^\Hi{t% (24) 
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where Koj{t' ; fi)=Tir[K{t' ] fi)pr]. Note that is now the interaction picture operator for the system due to 

the fact that [Hi{t'),Hr]=0. So far, all these derivations are exact. In order to attend to the complicated RHS of 
Eq. (j23ll . however, I resort to certain approximations. Following the standard concepts pertaining to the dynamics 
of open systems [S^, I make an important assumption, that Ki{t\ pl)k,Kqi{t\ for all itlj-terms on the RHS. 

This is justified if the field term, and the coupling term, V, are so very weak that the reservoir is almost 

not affected by these interactions. By imposing the conventional condition Tr[i?/(t')pr]=0 and performing a Markov 
approximation, I immediately obtain the following expression: 


rtj 


RHS= / dTSi{t')Si{T)Koi{t’;p){Ri{t')Rj{T))r 


ftf 


dTSl{t')Kol{t'-p)Slf^{T){Rl{t')Rl^{T))r 


- J dTSl{T)Kol{t'\p)Sl^{t'){Rl{T)Rl^{t'))r 
ftf 

+ J dTKoi{t';p)Sif^{T)Sifj,{t'){Rjfj_{T)Rif^{t'))r-, 


(25) 


where (• • • )r.=Trr[- • • Pr] are the reservoir correlation functions. If I further assume the RWA to be valid, the integrals 
in the above equation may be eliminated, and I arrive at a final form: 

RHS = -i[HLs,Koi{t';p)]-D*[Koiit';p)], (26) 

where H^s represents the Lamb shift [35j |. Considering that this is a standard procedure, I only highlight the key 
steps in Appendix II. If I transform Eqs. (IMD and (1261) back into Schrodinger’s theoretical framework and neglect the 
smaller Lamb shift, I reproduce Eq. ©■ 


The backward equation of the heat. As discussed earlier, the PDFs of the heat Qtem and Qqjt are the same in the 
master equation ©• Hence, defining one CF for the heat is adequate, e.g., ^{p). Previous results [i,[Tl[2i,|3il33 
have shown that this function can be evaluated by 

E{p) = = Tro[/3(t/; p)], (27) 

where p in the latter equation satisfies the same Eq. ([9]) but with a different initial condition of po. In comparison 
with the case of the work, one may expect that there exists an operator Eb(t';/i) (0 < t' < tf) that leads to the 
relation 


E{p)=Ti[Fo{0-,p)po]. (28) 

It is not difficult to find the backward equation of the operator if one follows either of the two methods related to the 
work. Here I only present the final result: 

dnFoit'-,p) = -Cfip)Foit';p) (29) 

with a terminal condition Fo(tf;p) = I. The current discussion also highlights the fact that, if the initial density 
matrix of the system is a completely random ensemble, pc, such as in the A^-level system, pc=I/N, the PDF of the 
heat must obey an exact equality, 

Ec[e"^'3QjT] ^ (30) 

where I used the subscript C to indicate the fact that all the QJTs star from this specific initial condition. The proof 
of this is reserved for Appendix HI. Note that there exists a classical version of the equality in the classical stochastic 
process 

I would like to assert that, since the use of the backward and forward equations related to the work (heat) always 
leads to the same CFs, from the viewpoint of computing, there are no specific reasons to favor one more than the other. 
This is only a question of personal taste or convention. On the other hand, when one carries out formal derivations, 
in some cases the forward method is more convenient than the backward method, or vice versa. For instance, one 
can easily obtain a concise expression related to the work equality if the backward method is used [H, [ 23 , HJl ■ 

Conclusion. In this paper, I have proved the equivalence of the two CF methods used for calculating work in weakly 
driven open quantum systems. Hence, in the present case, the PDF of the work defined using QJT is the same 
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as the PDF of the work defined for the combined system and reservoir using the TEM scheme. This finding has 
two implications. First, conceptually, the definition of the work and heat using the QJT is intuitive if one wants 
to understand these thermodynamic quantities from the viewpoint of the system. It is also not in conflict with the 
definitions that employ the TEM scheme on the system and reservoir. Second, since the recording of the QJTs has 
been realized in experiments [13, , it shall be realistic in practice to verify the work or heat equalities based on this 

notion. In addition, my current results could be extended to some specific quantum master equations. Eor instance, 
an immediate example is the adiabatically driven quantum master equations . Although concrete analyses are still 
to be conducted, the high similarity between this type of master equation and the current version strongly suggests 
this conjecture. Finally, I want to emphasize that conducting the TEM scheme on an arbitrary composite system is 
always possible, as this does not rely on the coupling strength or the time-evolution scales between the system and 
reservoir. However, whether or not QJTs exist for general cases, e.g., strongly coupling and nonMarkovian conditions, 
is still an unresolved issue. Hence, under these conditions, defining work or heat from the viewpoint of the system 
would become very challenging. I hope that new progresses in this direction can be presented in the near future. 


Acknowledgment. I thank Xin Wang for his participation in helpful discussions. I would also like to thank Editage 
(http:// WWW. editage.com) for English language editing. The work was supported by the National Science Eoundations 
of China under Grant No. III74025 and No. 11575016. 


APPENDIX I: AN EXPLANATION OF EQ. (IT^ 


I discuss here the time-dependent Hamiltonian of a closed quantum system, H{t)=Hs+Hi(t). This Hamiltonian, 
Hs, is assumed to have discrete eigenstates and eigenvalues: Hs\n)=€n\n). According to the TEM scheme [l^, one may 
define the exclusive work as WTEM=^n—^m, whereci, i=m, n denotes the energy eigenvalues of Hg that are measured 
at the beginning and the end of the nonequilibrium process. The CF of the PDF of the work is $(/i)=Tr(/i)ps], 
where 


Kt{p) = (31) 

U{t) is the time-evolution operator of H(t), and ps is the thermal density matrix of the bare system, Hg. One way 
to calculate $(/i) is to directly solve U(t). An alternative method involves determining an evolution equation about 
Kt{p). However, there is no such closed equation about Kt{p) with respect to the time, t. This problem may be 
circumvented by introducing 


= U{t')U'f{tf)e^f^^‘Uitf)U\t')e-^^^\ (32) 

It is obvious here that K{0; p)=Ktj{pL) and K{tf, p)=I. Interestingly, this new operator satisfies the following closed 
evolution equation about the backward time, t', 

If the bare system is composed of the aforementioned combined system and reservoir, and the interaction term, V, is 
assumed to be negligible, I arrive at Eq. m- 


APPENDIX II: SEVERAL KEY FORMULAS USED IN DERIVING EQ. (1^ 


The decomposition of the operator S implies that 


LJ UJ 


(34) 


Substituting these into Eq. (1251) and performing RWA, I obtain 


RHS = ^S{uj)S\uj)Koi{t'-n) / dse^‘^^Ri{0)Ri{-s))r 

u: “'O 

- V S{uj)Koiit'-' dse*‘"^(i?/(0)i?/^(s)), 

u: “'O 

- V S{uj)Koi{t'; f ' dse-^‘^‘^{Ri{0)Ri^{-s)) 

■’0 

+ y^Koiit';fi)Siuj)S^iu;) [ ' dse-“*(R/^(0)R/^(-s)),. (35) 

u: “'O 

If the correlation functions decay very quickly, these integrals can be approximated using one-sided Fourier transforms 
by replacing the upper limit tf—t' with infinity. Furthermore, it is useful to rewrite these one-sided Fourier transforms 
using the positive double-sided Fourier transforms 7 ( 0 ;) = {Ri{0)Ri{—s))r [^ : 

^“dse-(i?z(0)i?7(-s)). = ^jico) + (36) 

dse^^%Ri{0)Ri^{-s))r = (37) 

where V denotes the Cauchy principal value of the integral. After a simple rearrangement, I obtain Eq. ([26l) . where 
the Lamb shift term 


/ +00 

7(fl)/(a; - fl)dfl. 

-00 


(38) 


APPENDIX III: PROOF OF EQ. (1501) 

Two key observations are used here. First, Eqs. (EZl) and (1281) are in fact valid for any initial density matrix that 
is diagonalized in terms of the energy basis ®-g '5 

= Tr[Fo(0;/r)pc] = Tro[p(t/;/x)]. (39) 

The modified reduced density matrix p(t;/i) still satisfies Eq. ([9]) but with the initial density matrix of pc- Second, 
by applying the detailed balance condition about the rates and S{—oj)=S\uj)^ one may find 

A'(^/3)(pc) = 0. (40) 

With these two facts, I can easily prove the equality by using the backward equation (1^ or the forward equation ([^. 
For the former method, I need to evaluate the following expression: 

at/Tr[Fo(f';/x)pc] = -Tr[£*,(^)(Fo(t';Ai))pc] 

= -Tv[Fo{t'-p)Ct,{p){pc)] (41) 

Clearly, by substituting p=il3 into the above equation, performing an integral on both sides from time 0 to tf, and 
applying the terminal condition Foitf, p)=I, I immediately obtain the identity, 

Tr[Fo(0;x/3)pc] = l. (42) 

According to Eq. (IMl) . I may then prove Eq. (1501) . However, it is also easy to arrive at this equality if one applies the 
forward equation (0). I first re-express Eq. dSH) as 

Tro[p(t/;/x)] = TT[G(tf,0; p)[pc]], 


(43) 
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where the superpropagator G(t/, 0; fj.) equals T_ exp Cr{p)dT . Choosing and examining Eqs. (15^ and (HOI) 

again, I re-obtain the equality (p0|) . The reader is reminded that these discussions are also able to verify the validity 
of the work equalities [ 2 I, 
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